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Abstract 

It is known that in loop quantum cosmology (LQC) the universe avoids the singularity by a 
bounce when the matter density approaches the critical density pc (the order of Planck density). 
After incorporating the inverse volume modifications both in the gravitational and matter part 
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\^ • in the improved framework of LQC, we find that the inverse volume modification can decrease 

(N 

the bouncing energy scale, and the presence of nonsingular bounce is generic. For the backward 
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I evolution in the expanding branch, in terms of different initial states the evolution trajectories 
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, The other class with smaller initial energy density evolves back into the region a < a^h- In this 



classify into two classes. One class with larger initial energy density leads to the occurrence of 
bounce in the region a > a^h where a^h marks the different inverse volume modification region. 



region, both the energy density for the scalar field and the bouncing energy scale decrease with 



t^-' the backward evolution. However, in the deep modification region, because of the inverse volume 

op. 

modification the scalar field is frozen, such that the bounce is present when the bouncing energy 
• 1— I . 

^ ' scale decreases to be equal to the energy density of the scalar field. Using numerical method, we 



show the evolution picture for the second class bounce. 
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I. INTRODUCTION 



In cosmology, an outstanding problem is the big bang singularity which is expected to be 
solved by quantum gravity. Loop quantum cosmology (LQC) uses the framework developed 
from loop quantum gravity (LQG) 1] to deal with the issues in cosmology I2I, As a 
ramification of LQG, LQC inherits the nonperturbative and background independent feature 
of LQG. In LQC, the underlying geometry is discrete as in LQG, and the universe evolution 
is described by the difference equation which can go through the singularity nonsingularly 

The central technical issue in LQC is to quantize the Hamiltonian constraint which gener- 
ates the dynamical law. The classical Hamiltonian constraint consists of the basic variables, 
i.e., the SU(2) valued holonomies and the densitized triads. Conventionally, on quantization 
of the Hamiltonian constraint, the SU(2) valued holonomies in the gravitational part traces 
over the J = 1/2 fundamental representation, while in the matter part the SU(2) repre- 
sentation can be freely specified to define the inverse volume resulting in the quantization 
ambiguity known as inverse scale factor modification. Based on such quantum Hamiltonian 
constraint, the effective dynamics can be obtained by some approximation. By use of the 
effective dynamics one can investigate some phenomena in the semiclassical region, where 
the spacetime recovers the continuum and the difference equation is replaced by a differential 
equation A number of results, such as a natural inflation from quantum geometry 

7, ^, avoidance of a big crunch in closed cosmology 9|, an oscillating universe with suitable 
initial condition for inflation appearance of a cyclic universe 11 1 and a mass threshold 
of black hole [l^, etc, are interesting and remarkable. 

As shown in Ref. 13|, the quantization ambiguity also can appear in the gravitational 



part. Using arbitrary J representation for the holonomies, the Hamiltonian constraint op- 
erator is constructed in Ref. [l^ . It shows that in the semiclassical region the gravitational 
Hamiltonian gets modification similar to the inverse volume modification in the matter 
part. This is because that, in the classical Hamiltonian constraint, both the gravitational 
part and matter part contain inverse volume terms, and, for quantization, promoting the 
inverse volume into Poisson bracket between the volume and the holonomies leads to the 
same quantization ambiguity [14]. Therefore, as a phenomenal investigation, in the semi- 
classical region the inverse volume modification appearing in the gravitational part should 
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also be included. 



Recently, the semiclassical state is constructed in Refs. [15|, |16|, and the most important 
result is that, by evolving the semiclassical state backwards in the high energy density 
regime, the expanding universe bounces into a contracting branch, so that the singularity 



3; the discrete quantum 
16l . llTj l . This constraint 



can be avoided. An effective Hamiltonian constraint incorporatin 
geometry can well describe the evolution of the semiclassical state 
predicts a quadratic density correction in the modified Friedmann equation for the Hubble 
rate H = a/a, a p(l — p/ pc). The modified Friedmann equation implies a bounce when 
the matter density approaches the critical density pc (Pc is about 0.82 times the Planck 
density). Using the modified Friedmann equation, some interesting results are obtained: 
a bounce can happen avoiding of the singularity when the energy density approaches the 



critical value pc 18|]; the scaling solutions of the modified Friedmann equation have dual 
relationship with those in Randall- Sundrum cosmology [l^; the future singularity can be 
avoided by the modified Friedmann equation 2J]. However, in these works the inverse 
volume modification is neglected both in the gravitational part and the matter part, and 
for the gravitational part the holonomies are valued only on the J = 1/2 fundamental 
representation. Furthermore, from the viewpoint of lattice model, a large spin J means a 
longer range interaction between lattices which appear as background structures [20], and 



2l|. 



the inverse volume modification also could contribute to the seed of structure formation 
The effect of the inverse volume modifications to the effective dynamics is still unclear for 
the improved framework. Some further questions still need to be answered: What is the 
effect of the inverse volume modification in the gravitational part of the effective dynamics? 
What is its effect on the matter part? We shall look into these questions in the framework 
of the effective Hamiltonian. 



In Ref. 



IGj, the improved Hamiltonian constraint operator is introduced. The minimal 



area gap can be better exploited by the improved Hamiltonian to realize the physical idea, 
i.e., in the fiat model a bouncing universe is present purely because of the discrete quantum 
geometry. In this paper, we work in the improved framework to analyze the effect of the 
inverse volume modification both in the gravitational part and matter part of the effective 
dynamics. We find that the inverse volume modification in the gravity sector raised by a 
bigger quantization ambiguity parameter J could decrease the energy scale at the bounce 
point. And, whether a bounce really happens or not depends on whether the matter density 
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at the high energy region attains the bouncing energy scale. The inverse volume modification 
greatly changes the matter density below the scale a^. With the backward evolution in an 
expanding universe, in the region a < a^h the scalar field is frozen such that a nonsingular 
bounce is present. Our numerical results show the genericness of bounce in LQC. 

This article is organized as follows. In Sec. [ITl we give a brief review over the Hamiltonian 
constraint and the modified Friedmann equation. Then in Sec. Illlt the effect of the inverse 
volume modification in universe evolution is analyzed in detail. Section IIVI contains our 
numerical results, displaying the evolution picture of the nonsingular bouncing universe. 
Finally, the concluding remarks are made in Sec. |Vl 

II. HAMILTONIAN CONSTRAINT AND MODIFIED FRIEDMANN EQUATION 
A. Classical and quantum Hamiltonian constraint 

LQC is a canonical quantization of cosmology model mimicing the quantizing proce- 
dures as in the full theory. Imposing on the homogeneous and isotropic symmetry in LQC, 
the classical phase space is reduced to two degrees of freedom consisting of the conjugate 
connection c and triad p. The Poisson bracket between the conjugate variables satisfies 
{c, p} = |7/t, where k = SttG {G is the gravitational constant), and 7 is the dimensionless 
Barbero-Immirzi parameter whose value, 7 ~ 0.2375, set by the black hole entropy calcu- 
lation. For the fiat model, the usual Friedman-Robertson- Walker (FRW) metric expression 
can be identified with the relation given by 



where a is the FRW scale factor, and the absolute value of p denotes the two orientations of 
the triad. Here, we only take the positive orientation. In terms of the connection and triad, 
the classical Hamiltonian constraint is given by 



On quantization in LQC, there is no operator directly corresponding to the connection 
c. Instead, the elementary variables are the triad p and the holonomies hi{^) defined as the 
connection c along an edge, i.e., hi^n) = where /i is the length of the ith edge, and 



c = 



7d, IpI = a^ 



(1) 
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is a basis in the Lie algebra su(2) satisfying [xj, Tj] = EijuT^ . As in LQG, the holonomies and 
triads have well-defined operators. So, the Hamiltonian constraint given by Eq. ([2]) must 
be reformulated in terms of the holonomies and the triads. The gravitational part of the 
Hamiltonian constraint Hq in the full theory can be written as 

where is the curvature component of the connection. 

The strategies for quantization are: first, to reformulate the Hamiltonian constraint by the 
holonomies and triads, and then, to promote these elementary variables to the corresponding 
operators. The term E!^Ej / a/ | det E \ contains the inverse volume, and on quantization the 
inverse volume is promoted to the commutator of the holonomy and volume, resulting in the 
quantization ambiguity labeled by the spin representation J. For the J = 1/2 fundamental 
representation, the quantum operator for the gravitational part is given by 



^°=44^^''''^'' '''''' 

ijk 

6i . 2 ffic\ 2 ff^c 
— ■ sm — cos ' 



x(-(f)--(f)— (f)--(f))' 

where ip = Gh {h is the reduced Planck constant) and the volume operator V = p^^'^ . In 
Eq. (j4]) the holonomies are valued along the edges of the square loop whose area is fixed 
by the minimal eigenvalue of area operator in LQG. The physical area of the square loop is 
JjL^p = a£p, where a is order one. The area of the loop fixed by the minimal area of LQG 
shows the discrete feature of quantum geometry. We shall see that this point essentially 
captures the effect of quantum geometry making the effective dynamics very different from 
the classical case in the high energy regime. As for the arbitrary J representation, the 



gravitational Hamiltonian operator can be constructed as done in Ref. [13|. One only 
replaces the fixed parameter /io with p,. As shown in Ref. [igI, replacing /xq with /2 is a 
key step to realize the improved Hamiltonian constraint. After obtaining the gravitational 
Hamiltonian constraint operator for arbitrary J representation, it is useful to extract the 
effective theory which is expected to capture the key feature of the Hamiltonian constraint 
operator in the semi classical region. Next, we will show the inverse volume modification 
raised by the quantization ambiguity parameter J in the effective theory. 



B. Effective Hamiltonian constraint 



It is shown in Ref. [l^ that the quantum evolution determined by the quantum Hamil- 
tonian constraint is nonsingular in big bang point. By evolving the semiclassical states 
backwards, the universe evolution indicates that the big bang is replaced by a bounce when 
the matter density approaches the critical value pc- This feature can be well described by 
the effective Hamiltonian constraint in the region which is above the Planck scale and where 
spacetime recovers the continuum. The effective Hamiltonian can be obtained by the WKB 



approximation or by a systematic way 

arbitrary J representation is given by Q 



Here, the effective Hamiltonian for the 



H, 



^Sj{p) sin (pc) + H 



M: 



(5) 



where Sj (p) encodes the inverse volume modification for the gravitational part and is given 
by 

(6) 



Sj (p) = ^/pS { — 



and 



-I- 

Vg\lO L 



(g + 1) 2 + sgn (g - 1) |g - 1| 



(7) 
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\q-l 



where pj = ^^^y^£p mark the modification region, below which (i.e., q < 1) the modification 
is dramatic, and above which the effect of the inverse volume modification can be neglected 
with S (q) ^ 1. 



C. Modified Priedmann equation 



1. Small J modification 



One can assume that J is small, then Sj {p) ~ ^/p and the effective Hamiltonian becomes 

3 



H, 



^/p sin (pc) + H 



M- 



(8) 



Such an effective Hamiltonian presents a singularity-free bouncing universe with the modified 



Friedmann equation 



16 



18| 



3 V pc 



(9) 



where H is the Hubble parameter, the critical density pc = 3/^7^/2^0^ = and the 

matter density p = Hm/p^^'^. Based on the effective equation ([9]) some interesting issues 
have been discussed [18|, [l9|, \24\. 



2. Large J modification 

Now, we discuss the large J modification to the effective dynamics. By the effective 
Hamiltonian constraint given by Eq. the Hamiltonian equation can be obtained as 

dHeff 



P = {P, Heff] 

The Hamiltonian equation (fTOj) gives out 

1 



3 dc 



aa = -^Sj (a^) sin (pc) cos (pc) . 
7/i 



(10) 



(11) 



Squaring the above equation and then making use of the Hamiltonian constraint -f^e// ~ 0, 
the modified Friedmann equation incorporating the inverse volume modification can be 
obtained as 



P 



^P S{q)-^ 



Pc 



(12) 



where q = p/pj = {a/ag^Y 



III. EFFECT OF THE INVERSE VOLUME MODIFICATION IN UNIVERSE 
EVOLUTION 

A. Modification in the gravitational sector 



The modified Friedmann Eq. fll2p gives out the evolution of the universe. In the early 
time of the universe, if the matter density p = S {q) pc, then there exists a turn-around. 
At the turn-around, a bounce or a collapse happens depending on the sign of the second 
derivative of the triad with time. For the effective theory one can analyze the bounce 



behavior of the universe as discussed in Ref. 
c = {c, Heff}, one can get 

dlnSj {p) 



18| . From the Hamiltonian Eq. (|TOl) and 



P\p=0 



KSjip) 
nSjip) 



dlnp 
dlnSj (p) 
dlnp 



pc - 

Pc + 



2 1 OH. 



M 



S{q) \3^/P dp 
1 



P 



S{q) J 



(13) 




as in Ref. 25|]. For a 



p\p=0 = nS'^j{p)Pc o 1 + 



din Sj ip) 



dlnp 



) 



(14) 



For p < pj, the function d In Sj (p) /dlnp monotonically decreases, and when a -C a* = y/pj, 



that the recollapse requires < — 1, which classically violates the null energy condition and 
implies an expanding universe with increasing energy density. This exotic matter behaves 
like phantom field. In the next section, we shall explicitly give out the expression of the 
effective state parameter equation ujeff] it will show that the presence of a turn- around 
always leads to a bouncing universe. So, for the universe with quantum modification, the 
state parameter equation satisfies > — 1. 

In order to better understand the effect of the inverse volume modification on the modified 
Friedmann, we shall distinguish the two basic scales: One is the bounce scale abounce, which 
is determined by the condition p{a) = S (q) pc, and whose explicit value is related with the 
mater content (in the next section we will deal with the related issue). The other one is the 
inverse volume modification scale a* = ^/pj = ^JSn'^Jji/Mp = (8710^/^7/3)^''^ J^^^ip {j2 is 
replaced by the fixed area relation /i^a^ = a^p), whose magnitude depends on the parameter 
J. We know that abounce is the minimal scale for the evolution of the universe. If the inverse 
volume modification scale a* < abounce, the effect of the inverse volume modification should 
be neglected because in this range the inverse volume modification tends to its classical form 

(p) = y/P- However, for a large enough J, such that abounce < a < a^:, then the inverse 
volume modification becomes notable in LQC 

From the above analysis, we know that for the bouncing universe the matter density at 
the bouncing energy scale is Pbounce = S (g) pc- In the region abounce < a < a^, the function 
S (g) monotonically increases and S {q) < 1. For abounce <^ a*, 5 (g) ~ |^ <^ 1. Therefore, 
we can conclude that for LQC the inverse volume modification in the gravitational part 
reduce the matter density at the bounce point. 

In this subsection, we have analyzed the effect of the inverse volume modification to the 
bouncing universe. However, for LQC, whether a bounce happens or not depends on the 
content of the matter field. In the next section, we will show the effect of the inverse volume 
modification in the matter sector on the evolution of LQC. 



P\v=o 



K.Sj{p)pc [^ + uj); a a*, p\p=o npjpc {1.05 + u). Therefore, one can obtain 
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B. Modification in the matter sector 



In this section, we focus on the inverse volume modification in the matter sector for 
the modified Friedmann equation. In earlier literatures, the SU(2) representation for the 
matter part is labeled by a half integer which can be freely specified independent of the J 
representation in the gravitational part. As argued in Ref. 13|], it is more natural to take 
the same representation for them. Therefore, for a massive scalar field with self-interaction 
potential V (0), the Hamiltonian Hm is given by [tI] 

1 



M 



-dj (a)pl + aV(0), 



(15) 



where is the conjugate momentum for the scalar field 0, and dj (a) is the eigenvalue of 
the inverse volume operator which can be approximately described as 

.2- 



dj (a) 



D{ — \a 



(16) 



In which, the modification region is still marked by the scale a, = (Svra^/^^/S)^^^ J^/^£p, 
and the function D {o? /al) respects the inverse volume modification and has the form 



D{q) = [j^] g^|7[(g + l)^-|g-l| 



llg 



7 ^1 1 ^ 

{q+l)'^ — sgn{q—l)\q — l\'^\> . 



(17) 



Above the scale a*, D (g) ~ 1, so the inverse volume modification can be neglected. While 
at g < 1, D (aVa') ~ (12/7)^ (a/a*)'^ 

Now, the modified Friedmann equation becomes 



K 



Pc\ 



(18) 



where the matter density p = Hm/p^^"^ which receives the inverse volume modification. 
From the Hamiltonian equation = {0, Hm}-, one can get 



p^ = d/ (a) 0. 



(19) 



It leads to the matter density 



+^(0) 



(20) 
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Differentiating Eq. (fT9i) and substituting it into another Hamiltonian equation = 
{p^,Hm}, the modified Klein-Gordon equation can be obtained as 

i=-3n(l-~)4'-DV(^), (21) 



3 (i In a 

where prime denotes Then the time derivative of the matter density is 

-3^-^644^1 (22) 



6 dlna 

where din D/ dlna is a monotonically decreasing function in the region a < a*. So, there 
exists a characteristic scale ach ~ 0.91a* determined by din D/ din a = 6, and the mater 
density attains the maximal value. Classically, with the expansion of the universe, the 
matter density decreases. However, in the region of a < ach, the matter density increases; 
above this region the matter density starts to decrease and enters the regime where the 
quadratic density modification dominates and the effect of the inverse volume modification 
can be neglected. In the following, we will discuss the possible evolution of LQC in different 
regions. 

First, it is useful to identify the effective matter density and pressure. From Eq. f|T8l) . 
the modified Raychoudhuri equation is written as 

d 



H + H 

a 



2 



Then the effective energy density and pressure can be defined as 



p. 



Peff = p(S~^], (24) 



They satisfy the conservative equation 

Pcff + 3H (pe// + Peff) = 0. (26) 

Substituting Eq. ( l22l) into the expression of Peff, the effective state parameter equation 

is 

dS 



peff 02 + 2DV {(p) ^- ~ ^ 
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where the function dS/dq monotonically decreases in the region a < a*, and asymptotically 
vanishes beyond this region. From the effective parameter equation we know that both 
the inverse volume modification and the quadratic density correction can contribute to the 
inflation of the universe. 



26i. So, 



In LQC, the evolution trajectory is determined by the dynamical initial condition 
we evolve the universe at large scale a ^ a* backwards in the expanding universe where the 
matter density p <^ pc and the quantum modifications can be neglected. With the backward 
evolution, the matter density increases in the region a > a^h- If the matter density satisfies 
p ^ S (alh/al) Pc then pe// = leads to a bouncing universe which is caused purely by the 
quadratic density modification. This situation is similar to the case of dftomce ^ O'*, where 
the inverse volume modification can be neglected. If the mater density p < S{a1f^/al)pc, 
the universe evolves back to the deep inverse volume modification region a < ach- In the 
region a < ach, from Eqs. ([7j) and fl22|) we know that both the matter density p and the 
bouncing energy scale S{q)pc decrease with the backward evolution. If the bouncing energy 
scale S{q)pc decreases faster than the matter density, then there is no bounce; the universe 
will evolve towards the big bang point where the effective theory is invalid and the quantum 
difference equation takes over the evolution of the universe. Another possible way is that 
the mater density decreases more slowly than the bouncing energy scale, until p = S{q)pc, 
where a bounce happens and the universe is bounced into a contracting branch. We would 
like to point out that the actual evolution way favors the latter. It is because that, with the 
backward evolution, the density decrease. However, for the density p = |-D^V^ + ^ (0); in 
the small volume, the term D^^ cx (a/a*)^^ is far larger than one. This needs the velocity 
to be extremely small, such that the field is frozen. So, the kinetic term can be neglected, and 
the potential dominates. With the decreasing scale a, the field is frozen, until the bouncing 
energy scale S{a?/a^)pc decreases to be equal to the potential V^(0). At that point, a 
bounce happens. In what follows, using the numerical way, we will show the evolution of 
the universe. 

IV. GENERICNESS OF BOUNCE SHOWN IN NUMERICAL RESULTS 

In the inverse volume modification region, a negative scalar potential can lead to a bounc- 



ing universe only when the condition p = is satisfied ll|, |27|. This is also true for the 
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modified Friedmann equation presented in this paper. From the modified Raychoudhuri Eq. 
fl23|) . it is easily seen that in the region a < a^h the vanishing density imphes a bounce. Such 
bounce purely comes from the inverse volume modification. The reason is that the inverse 
volume modification in this region makes 



a, 1^ T f IdhiD 
- d=o = -77^ V 1 - 77-71 



which implies a bounce. In the following, we will show the evolution trajectory of LQC with 
a positive potential. 

Using the modified Friedmann Eq. ( ITSl) and the Klein-Gordon Eq. ( 1211) . we can obtain the 
dynamical law for the universe with a massive scalar field with the potential V (0) = m^0^/2. 
Here, for the numerical calculation we take the mass m = 0.2. From the above discussion 
we know that these trajectories can be divided into two classes: the ones with higher energy 
density can be bounced into an contracting branch in the region a > a^h, and the others 
with lower energy density p < S {alh/al) pc ~ 0.94pc evolve back to the deep modification 
region a < ach- In this paper, we are only concerned about the latter. For the backward 
evolution, we take the initial state marked by the energy density at the scale ach- Our 
numerical results show that the presence of a bounce is generic for theses trajectories with 
the initial state pm in the region (0.001 ~ 1) 5* {alhl'^D Pc- Below this region, we think that 
it is not suitable to take the energy density far below the Planck density as the initial state 
in the quantum modification region. The initial state can be written as 

pin (ach) = ^D'^ {al,^/al) (p"^ + ^m^^^ 

= XS [alJal] Pc, (28) 

where A denotes the initial state with different initial energy. From Eqs. (fT8|) and (pTj) 
we know that under the rescaled transformation a a/a^ the dynamical equations do not 
change and do not depend on the scale a explicitly. So, under this rescaling we can take an 
arbitrary value of J to show the evolution trajectories. For other choices of J, the evolution 
trajectories are the same as the chosen one. Here, we take J = 1000 to show the evolution 
picture. 

Figure [1] shows the presence of bounce for A = 0.5. For different initial values of (j), 
the bouncing scale is different. With larger value of |0|, the bouncing scale is larger. The 



contrary is that with larger value of 



the bounce happens at deeper region. This can 
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be regarded as that, for the backward evolution, a larger kinetic term takes longer distance 
to be frozen. At the bouncing point, the field is frozen at some fixed value of (j). Figure [2] 
shows that with the backward evolution the scalar field is frozen and then bounced into a 
contracting branch. For other choices of A, the bouncing behavior is similar as presented in 
Figs. [U and [H The difference is that with smaller value of A the bounce happens at deeper 
region. 

Furthermore, in the inverse volume modification region the scalar field can be pushed 
to climb up the potential hill. As a result, this can set the initial condition for a classical 
inflation 28|. To achieve 60 e-folds of inflation for a quadratic potential V = m^0^/2 , the 
field must be displaced by at least 3Mp; from the minimum of its potential at the onset of 
inflation. For the case of A = 0.5 presented in the above, this means that for a successful 
inflation the bounce should happen at bigger scale for larger value of \(f)in\ shown in Fig. [H 

V. CONCLUSION 

In this paper, we mainly discuss the effect of the inverse volume modifications on the 
effective dynamics, which implies a nonsingular bounce for the evolution of the universe 



when the matter density approaches the high energy region. This is a nontrivial resu 



compared with the bouncing mechanism in the fundamental J = 1/2 representation Id Il8|. 



t 



Because until now there is still lack of a complete quantization for large J representation 

ution picture presented in this paper 



16| . Furthermore, we argue that the 



in the improved framework, we expect that the evo 
will be given by a complete quantum theory as in 
genericness of bounce in the high spin representation is not accidental: the presence of 
bounce is a consequence of the fact that LQC essentially incorporates two quantum effects, 
i.e., the inverse volume modifications and the quadratic density modification raised by the 
minimal area gap, and the two quantum effects jointly induce the occurrence of bounce. 
For the inverse volume modification, the modification region is marked by the scale a^,, 
below which the modifications are notable and above which the modifications become weak. 
The inverse volume modification in the gravity sector decreases the bouncing energy scale 
in the region abounce < a < cig* (neglecting the inverse volume modification the bouncing 
energy scale is equal to the critical density pc )■ And, the inverse volume modification in 
the matter sector changes the classical behavior of density. Below the scale ach ~ 0.91a* the 
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FIG. 1: The evolution of the scale factor is shown for a scalar field with potential V {(f)) = m?4P' /2 
for different initial states distinguished by different initial values of (these initial states are 
constrained by the fixed initial energy density = 0.55 (cl\/'^'*) P'^ ~ O.YApc at the scale a^h, and 
(pin is determined by Eq l28p . The solid lines denote the expanding branch, and the dashed lines 
show the contracting branch. In the region a < act ~ 0.91a*, the backward evolution shows that 
for bigger initial value of {(pml the bounce happens at bigger scale. 
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FIG. 2: The sketch shows that near the bouncing scale the scalar field is frozen such that the 
kinetic term vanishes and the potential dominates. Here, the solid lines and the dashed lines 
denote, respectively, the expanding branch and the contracting branch. The initial states take 
values as in Fig. 1 constrained by the fixed initial energy density pi„ = 0.5S (cic/i/^*) Pc ^ 0.74/9c- 

time derivative of the matter density is positive, implying that, with the backward evolution 
of an expanding universe, the matter density decreases. So, in the region a < ach the 
presence of bounce depends on the relative decreasing rate between the matter density and 
the bouncing energy scale. However, in the deep quantum modification region the scalar 
field is frozen. So the kinetic term can be neglected, and the matter density is dominated 
by the potential. Therefore, the matter density is almost frozen at a fixed value. With the 
backward evolution, a bounce happens until the bouncing energy scale decreases to be equal 
to the potential. Our numerical results shown in Figs. [1] and [2] provide a illustration of the 
evolution pictures. For different initial values of the matter density the bouncing scale is 
different. 

In earlier literatures, the SU(2) representation for the gravitational Hamiltonian is taken 
to be the fundamental representation, but for the matter Hamiltonian the representation is 
freely specified. If one does so, then in this paper the modified Friedmann equation becomes 
= /tp (1 — p/ Pc) /3, where p = \D~^(j? + V (0) receives the inverse volume modification. 
For the backward evolution in the expanding branch, if the matter density increases to be 
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equal to the critical density pc in the region a > ach, a bounce happens. If the universe 
evolves back through the scale ach, the matter density starts to decrease, so a bounce will 
never be present for the universe with a positive potential. This is a distinct difference 
between the two views of the inverse volume modification. Moreover, the inverse volume 



modification can leave indirect imprint on the CMB spectrumj28|. We expect the further 
research could put constraint for the value of J based on the prediction for observations. 

In this paper we analyze the effect of the inverse volume modifications on the improved 
LQC, and we show that the presence of nonsingular bounce is generic. 
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